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O . Abstract 

Q ■ We study rapidity gaps between jets in pp collisions at the Tevatron by a novel 

, solution of the nonforward BFKL equation including nonleading effects through the 
consistency constraint and running coupling. Results differ from earlier calculations 

pi • based on the asymptotic Mueller- Tang formula, but agree with DO data when includ- 

D . ing full event simulation using PYTHIA to model the gap survival probability. 
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Mueller and Tang [||] (MT) proposed that events with two forward-backward high-i^T jets 
with a large rapidity gap (of size y) between them could provide a test of the nonforward 
BFKL equation, in the form of elastic parton-parton scattering by hard (|t| ^ 0) colour 
singlet exchange. However, the predictions obtained for the t and y dependence are not 
confirmed by D0 and CDF data. To crudely reproduce the t-shape of the D0 data 
the running of as{—t) had to be disregarded and the y dependence measured by CDF 
remains unexplained. The MT approximation is valid only for asymptotically large ener- 
gies and neglects some contributions which may be relevant for the present experimental 
conditions. Here, we show that a more complete solution of the BFKL equation gives a 
good description of the data. 

The cross section da/dt for elastic parton-parton scattering via BFKL pomeron ex- 
change is dominated by the imaginary part of the amplitude 

lmA(v t) = I ^o(k,Q)$(^,k,Q) 

^ J vr [(k + Q/2)2 + so][(k-Q/2)2 + 5o]' ^ ^ 



^Talk presented by RE at the 9th International Workshop on Deep Inelastic Scattering, Bologna, April 
2001, to appear in the proceedings. E-mail: enberg@tsl.uu.se. 
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Figure 1: Parton level gap fractions (relative to LO QCD one-gluon exchange) in quark- 
quark scattering from our BFKL solution and the asymptotic MT approximation with the 
running coupling as{—t) in the prefactor, versus gap size y (for the given E'^- values) and 
versus Et (for the given values). In the right-hand plot, the MT curve is multiplied by 
a factor of 20. 

where Q/2 ± k denote the transverse momenta of the exchanged gluons and t = —Q^. 
The scale sq is a physical infrared cut-off that reflects the confining properties of the QCD 
vacuum which suppresses the propagation of gluons with small virtuahties. We set sq to 
1 GeV^, but note that the particular value affects mostly the normalisation and not the y 
and t dependence. 

The impact factor $o(k, Q) is an effective vertex between the pomeron, the proton 
and the final state containing the high Ex parton(s). At large \t\ and y, the pomeron 
predominantly couples to individual pointlike partons |Q, corresponding to 



We assume that this is a good approximation also in the nonasymptotic region. $(?/, k, Q) 
is the gluon distribution including full y evolution which is found by solving the nonforward 
BFKL equation with $o(k, Q) as the input. The nonleading corrections to the BFKL 
kernel are taken into account in an approximate way by including the consistency constraint 
and using a running coupling where the scale /x^ = /c^ -|- is evaluated locally in the 
gluon ladder. 

The BFKL equation is solved numerically by expanding it in a truncated Fourier series 
and discretising the system by the Tchebyshev interpolation method [§. The resulting 
theoretical gap fraction is shown in fig. |I] together with the result of the MT formula which 
has a completely different behaviour. This is because MT corresponds to the asymptotically 
leading part of the scattering amplitude, which is not sensitive to the low scale details of 
the hadron structure and the QCD vacuum. However, the amplitude also contains a 
nonasymptotic part which depends on the physical cut-off scale and which dominates at 



$o(k,Q)~ v'a,((Q/2 + k)2 + So) a.((Q/2 - k)^ + sq). 
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Figure 2: Fraction of jet-gap-jet events over all jet events versus jet separation At] for 
15 < Et2 < 25 GeV and of the lowest energy jet; D0 data compared to our 
BFKL-based model implemented in pythia relative to all LO QCD 2 — 2 processes. 



small y. This part contains a log(Q^/so), visible at the Born level, which gets suppressed 
when y increases. This is due to gluon reggeisation 0, i.e., the negative virtual correction 
terms in the BFKL kernel. The numerical value of log((5^/so) is large and the BFKL 
evolution is not rapid enough to reach the high y asymptotic limit in the currently accessible 
domain of y. We stress, that this log(Q^/so) in not an infrared artifact, since it would 
also appear if we considered the analogous process of heavy onia scattering, where Sq 
i^QCD ^ So ^ Q"^) would be a hard scale corresponding to the heavy onium size. 

In order to compare with data, the parton level cross section discussed above must 
be convoluted with QCD evolved parton distributions of the incoming hadrons and the 
transition to final state hadrons be made. The latter is particularly important when con- 
sidering rapidity gaps, which are defined in terms of a region without particles. One must 
then consider all additional activity in an event, such as higher order parton emissions, 
multiple parton scattering and hadronisation. A simple way of doing this is to multiply 
the cross section with an estimated gap survival probability. However, this probability 
cannot be independent of the subprocess, with its associated colour string topology, or of 
the kinematics in the event. For example, gluon-gluon scattering and a larger hard scale 
resulting in more parton radiation is expected to give a smaller gap survival probability. 

The only way to handle these complex processes is via Monte Carlo event simula- 
tion. We have therefore implemented our cross section for elastic parton-parton scattering 
via colour singlet exchange as a new hard subprocess in the event generator pythia [Q. 
Thus, higher order parton emissions are included through conventional parton showers 
and the underlying event is modelled with multiple interactions (using the best MI model 
MSTP(82)=4). To be consistent with our description of rapidity gaps between a hard 
system {W, jets, bottom, J/ip) and a leading proton observed at the Tevatron, we also 



include the soft colour interaction (SCI) model |TT]]. The combination of SCI, with colour 



3 



exchange probability P = 0.25, and MI used here to model the underlying event is consis- 
tent with our previous study [0 . Although SCI was introduced to explain the creation of 
forward rapidity gaps, it has here the opposite effect by rearranging strings to span across 
the gap produced at the parton level by the BFKL singlet exchange. 

The gap fraction resulting from this complete model describes quite well the D0 data, 
as shown in fig. ^ The normalisation of the gap ratio is, however, uncertain for two 
reasons. The parameter sq affects the normalisation of the parton level cross section and the 
simulated survival probability depends on the SCI model parameter P. A more complete 



study, including CDF data, will be presented in a coming paper . 

In conclusion, we have solved the nonforward BFKL equation for colour singlet exchange 
including nonleading corrections. The part dominating at Tevatron conditions is not the 
asymptotically dominating part contained in the Mueller- Tang formula, but a part governed 
by a large logarithm that disappears in the large- ?/ limit. Using pythia with full event 
simulation we obtain a dynamically varying gap survival probability resulting in a good 
description of the Tevatron data on rapidity gaps between jets. 

We thank J. Kwiecihski for stimulating discussions and acknowledge the support from 
the Swedish Natural Science Research Council. LM is supported in part by the KEN grant 
number 5 P P03B 144 20. 
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